adp header will be provided by the publisher 



Superfluid to Mott-insulator transition in an anizotropic 
two-dimensional optical lattice 



**2 



T. P. Polak* 1 and T. K. Kopec 

1 Adam Mickiewicz University of Poznah, Faculty of Physics, Umultowska 85, 61-614 Poznan, Poland 

2 Institute for Low Temperatures and Structure Research, Polish Academy of Sciences, POB 1410, 50-950 
Wroclaw 2, Poland 

Received 22 October 2008 by U. Eckern 

Key words Optical lattice, Bose condensation, Mott-insulator transition. 
PACS 05.30Jp, 03.75.Lm, 03.75.Nt 

We study the superfluid to Mott-insulator transition of bosons in an optical anizotropic lattice by employing 
the Bose-Hubbard model living on a two-dimensional lattice with anizotropy parameter k. The compressible 
superfluid state and incompressible Mott-insulator (MI) lobes are efficiently described analytically, using 
the quantum U(l) rotor approach. The ground state phase diagram showing the evolution of the MI lobes 
is quantified for arbitrary values of k, corresponding to various kind of lattices: from square, through 
rectangular to almost one-dimensional. 



1 Introduction 

The type of order that a physical system can possess is utterly affected by its dimensionality. In two- 
dimensional (2D) systems with a continuous symmetry long-range order is destroyed by fluctuations at 
a finite temperature (Tj. However, the competition between ground states in T = can lead to a zero- 
temperature phase transition, driven solely by quantum mechanical fluctuations. In this context unconven- 
tional behavior in low-dimensional systems was intensively studied in the past years J2). 

Since the experimental realization of Bose-Einstein condensation many properties of it have been 
elucidated 015]. Atomic gases allow clean and controlled observation of variety physical phenomena 
from condensed matter physics, e.g., the Berezinskii-Kosterlitz-Thouless (BKT) phase transition with the 
emergence of topological order Q. The merging of atomic and condensed matter physics has opened 
exciting new perspectives for the creation of novel quantum states. Especially, systems of cold atoms in 
optical lattices |7, 8| facilitate an experimental environment, where a rich variety of quantum many -body 
models can be implemented in a wide range of spatial dimensions, geometries, and particle interactions. 
Among these topics the emergence of condensation and superfluid order in an optical lattice has been a 
major issue in recent years J9] [TOj QT) . In this context, the presence of the optical trapping structure offers 
a unique way to increase (or decrease) the dimensionality of the system in a clean experimental setup, thus 
providing a playground for studying the effect of dimensional crossover on the quantum phase transition. 
This can happen in optical lattices in which atoms can tunnel easily along one spatial direction but not 
along the other one. This motivates the analysis of very interesting physics, e.g., of an anisotropic array of 
coupled one-dimensional (ID) Bose gases. Here, the coupling is provided by the intersite tunneling that 
can be made variable by adjusting the optical lattice potential. A one-dimensional situation is created by 
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suppressing tunneling in two directions by using two standing waves with very high laser intensities that 
control the barrier between the lattice sites. 

The aim of this work is to study the superfluid (SF) to Mott-insulator (MI) transition by means of 
the Bose-Hubbard model in 2D optical lattices with variable lattice anizotropy parameter. In order to 
quantify the evolution of the ground state phase diagram, various kind of lattices, from square, through 
rectangular to almost one-dimensional will be considered. To this we have adopted a theoretical approach 
to the strongly interacting fermions lfl2l in the Bose-Hubbard model in a way to include particle number 
fluctuations effects and make the qualitative phase diagrams more quantitative iTPTl . The key point of our 
approach is to consider the representation of strongly interacting bosons as particles with attached "flux 
tubes". This introduces a U(l) phase variable, which acquires dynamic significance from the boson-boson 
interaction. In the present work we investigate the dimensional crossover from two- to almost ID lattices 
and obtain the ground state phase diagram. 



2 The model and method 

We start with the generic model for the Mott-insulator transition, namely the Bose-Hubbard model [ 14 1 



(1) 



where a\ and aj stand for the bosonic creation and annihilation operators that obey the canonical commu- 
tation relations [on, at] — Sij, rii — a\ai is the boson number operator on the site i, U > is the on-site 
repulsion and the chemical potential /i controls the number of bosons. Here, (i, j) identifies summation 
over the nearest-neighbor sites. Furthermore, tij is the hopping matrix element with dispersion 



t (k) = 2t (cos k x + k cos ky) , 



(2) 



where k is the anizotropy parameter and t sets the kinetic energy scale for bosons. By varying quantity k 
between zero (ID) and one (2D), different anizotropic rectangular lattices emerge. 
We write the partition function of the system 



[VaVa] exp — / drTi (r) — / drdi (r) — — a (r) 

, Jo v T 



(3) 



using the bosonic path-integral over the complex fields (r) depending on the "imaginary time" < r < 
(3=1/ k-oT with T being the temperature. We decouple the interaction term in Eq. (fl]i by a Gaussian 
integration over the auxiliary scalar potential fields 



Vi (r) = Via + V{ (t) 
with static [3Vm = Vi (lo v = 0) and periodic part 



+oo 



Vi (^~) = ft 1 J2V KO ex P (i^r) + c.c, 



(4) 



(5) 



where lo v = 2m>/(3 (v = 0, ±1, ±2, ...) is the Bose-Matsubara frequency. Periodic part V( (r) 
V( (r + (3) couples to the local particle number through the Josephson-like relation 



4>i (r) = V( (r) 



(6) 
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where <f>i (r) = dcfii (r) /dr. The quantity <fi (r) is the phase field satisfies the periodicity condition 
4>i (0) — 4>% (0) as a consequence of the periodic properties of the V( (r) field. Further, we perform 
the local gauge transformation to the new bosonic variables 

a% (t) = hi (r) exp [ifa (r)] (7) 

that removes the imaginary term — i drcfti (r) n, (r) from all the Fourier modes. From the above we 
deduce bosons have a composite nature made of bosonic part 6j (r) and attached "flux" exp [ifc (r)]. Note 
that a similar method was used in a functional-integral formulation to treat the quantum dynamics of a 
microscopic model of a Josephson junction, including the dissipative effects of quasiparticle tunneling 
fT31 . Next, we parameterize the boson fields 

bi (r) = & + b\ (r) (8) 

and incorporate fully our in calculations the phase fluctuations governed by the gauge U (1) group and 
drop corrections to the amplitude by assuming bi (t) = bo, which was proven to be justified in the large 
U/t limit we are interested in ifPJl [TH . By integrating out the auxiliary static field Vto we calculate the 
partition function Z = J [Dcf)] exp {— <S p h [<fi]} , with an effective action expressed in phase-only terms 



s ph [4>] = J^dr 



( r ) + l£0. ( T ) 



*W J.„e^*( T ) } . (9) 



We note that the phase action for the Bose-Hubbard model is closely related to the standard model of 
Josephson junction arrays, which contains the charging energy as well as the Josephson coupling Ifl6lll7| . 
The phase stiffness coefficient is given by Jy = bgiy , where the value 



1 / 1 



b '° = u (iv^> + /1 l (10) 

(id) 

is obtained from minimalization of the Hamiltonian &H (bo) jdbo = 0; Jl/U — n/U + 1/2 is the shifted 
reduced chemical potential. The total time derivative Berry phase imaginary term in Eq. ((9} is nonzero due 
to phase field configurations with 

& (/?) - <j>i (0) = 2nm z (11) 

where, nii — 0, ±1, ±2.... Therefore, we concentrate on closed paths in the imaginary time (0, (3) labelled 
by the integer winding numbers m^. The path-integral 

\v<j>}...= y j / [v<t>m / P4>{t))..., d2) 



n 



*((«) 



includes a summation over m, and in each topological sector the integration goes over the gauge potentials. 
To proceed, we replace the phase degrees of freedom by the unimodular scalar complex field tpi which sat- 
isfies the quantum periodic boundary condition tpi ((3) = tpi (0). This can be conveniently done using the 
Fadeev-Popov method with Dirac delta functional resolution of unity Ifl9l , where we take ipi as continuous 
but constrained (on the average) variable to have the unimodular value 

1 = J lVi>V$\ 5 (r)| 2 - N^j Y[5 (Vi - e^) 8 - e -**M) . (13) 
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Fig. 1 (online colour) The phase boundary between the Mott insulating and superfluid phases for the 
anizotropic lattice with k = 1, 0.5, 10 , 10 -2 , 10 -5 from left to right. Inset shows the tip of the first lobe 
(n B = l)forK = 10" 5 . 

Introducing the Lagrange multiplier A, which adds the quadratic terms (in the ipi fields) to the action Eq. 
(O, we can solve the constraint. Using such description is justified by the definition of the order parameter 

*s = (at (r)> = (h (r) exp [ifa (r)]) = b ^, B , (14) 

which non- vanishing value signals a bosonic condensation (we identify it as superfluid state). Note that 
a nonzero value of the amplitude bo is not sufficient for superfluidity. To achieve this, also the phase 
variables, must become stiff and coherent, which implies ipB 7^ 0. 

3 Phase diagram 

The partition function is written in the form 





~v\~ 


e -Nf3F{\) 


/ —too 







(15) 
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with the free energy density T = — In Zj j3N given by: 



^ In | [V^Vj] exp | ]T J dTdT ' K JI « 



+ X5ij) 5(t — t')] - jijipiipj 



where 2^ — 1 if i, j are the nearest neighbors and equals zero otherwise, 

Hj t ') = < cx p {-i y>i ( t ) - <t>j (t')]» 



(16) 



(17) 



is the two-point phase correlator associated with the order parameter field, where (...) is the averaging with 
respect to the action in Eq. ([9JI. The action with the topological contribution, after Fourier transform, we 
write as 



where 



(u v ) = A - J k + 7 1 (u v ) 



(18) 



(19) 



is the inverse of the propagator. The final form of the correlator, after Fourier transform, can be written as: 



1 4 V- 6XP 



.M 



m + 



1 - 4 (to 



it 

!7 



where 



* = £ 



exp 



-£//? (m + fl/U) /2 



(20) 



(21) 



is the partition function for the set of quantum rotors. The form of Eq. d20b assures the periodicity in the 
imaginary time with respect to fi/U + 1/2 which emphasizes the special role of its integer values (see 
Fig. [T]). Within the phase coherent state the order parameter ipB is evaluated in the thermodynamic limit 
N — ► oo by the saddle point method 8T '/ 6 A = takes the form: 



1 



iV/3 A - J k + 7-1 (w„) 



(22) 



The phase boundary is determined by the divergence of the order parameter susceptibility Fk=o (k>„=o) = 
0, which determines the critical value of the Lagrange parameter A = Ao that stays constant in the whole 
ordered phase. We introduce the density of states 



p«,«) = ^E*K-*( k )/ < ] 

k 

for rectangular anizotropic lattice (see inset Fig. |2j: 

i ... / l(i + K) 2 -e 



(23) 



4k 



[e (« - ie - 1|) + e (« - ie + 1 

20(1 -K- ' 



K 



4k 



7r 2 v / (1 + K) 2_ C 2 VV + 



(24) 
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Fig. 2 (online colour) The position of the tips of the lobes as a function of the anizotropy parameter k 
for fi/U — 0.5, 1.5, 2.5, 3.5, 4.5, 5.5, 6.5 from right to left in double logarithmic scale. Inset shows the 
evolution of the density of states p (£, k) from square lattice (k = 1) through rectangular (k = 0.33, 0.66) 
to one-dimensional (« = 0). 



where 9 (x) is the unit step function and K (x) is the elliptic function of the first kind [20|. With help of 
the above and after summation over u u , the superfluid state order parameter becomes 

l-rt-if /«■->« - (25, 

In Eq. dZSb u ([i/U) = frac (fi/U) — 1/2, where frac (x) = x — [x] is the fractional part of the number 
and [x] is the floor function which gives the greatest integer less than or equal to x. The zero-temperature 
phase diagram of the model can be calculated from Eq. (l25l >. We recover the results for pure 2D case 
(the lowest lobes in Fig. Q3 fl3l . Still there is a particle-hole asymmetry visible not in the position of the 
maximum of the lobe - like in simple cubic lattice - but in the shape of the curves. However, when we 
approach ID case the particle-hole symmetry is restored. In the context of the paper [21 J the existence 
of asymptotic restoration of the statistical particle-hole symmetry in ID dirty boson problem seems to be 
apparent. Next, we resort to the unimodular-field description and calculate the effects of the fixed boson 
number 

nB = n ^ (Sl (t) ai (r)) (26) 

i 

in the system. For the interacting problem, with the full phase action Eq. (0 we get 

_ . tie (A) within MI phase 

ns = <! n B (\ ) -2t/>%v(§) within SF phase ' ( } 
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Fig. 3 (online colour) Boson occupation number nj at T = for the anizotropic lattice with k — 10 
as a function of chemical potential fi/U and hopping t/U. The Mott insulator is found within each lobe of 
integer boson density. Inside each of MI lobe the integer occupation number nj is indicated. 



where ijjs is given by Eq. d25l ). In the limit T — ► an analytical solution of the total boson density consists 
of the occupation number for neutral bosons n& and a contribution Sni, from a fluctuating phase field: 

nn(X) = 77+2-HZ7 

, ^ (28) 

with SX = A — Ao- Here, the parameter A is self-consistently determined via Eq. d22l . Changes in 
the anizotropy k provide the condition for emerging MI (see Fig. [2]), with recognizable steps-like structure 
(Fig. [3). We cannot obtain pure ID case since even in T — long-range order is destroyed by the quantum 
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fluctuations. Regarding the comparision of our method with the previous appraoches, e.g., |[T4l[T8ll we note 
that the qualitative shape of the lobes resulting from our apprach is not the same for 2D and 3D cases, and 
steeper for the two-dimensional system. Furthermore, we found 1 1 3 1 that our results are in good agreement 
with the recently published quantum Monte Carlo calculations on three-dimensional Bose-Hubbard system 
El. 

4 Final remarks 

In conclusion, we have performed a study of the superfluid to Mott-insulator transition of bosons in an 
optical anizotropic lattice by employing the Bose-Hubbard model living on a two-dimensional lattice with 
anizotropy parameter k. The compressible superfluid state and incompressible Mott-insulating lobes are 
efficiently described analytically, using the quantum U(l) rotor approach. Our motivation comes from the 
fact that the development in experimental techniques of trapping and controlling quantum gases allow to 
investigate such phenomena by offering a platform for the exploration of highly non-trivial quantum phases 
and critical phenomena in dimensionality-tunable systems. The technique used in this paper can be easily 
extended to more general situations, including, e.g., multi-species bosonic systems. Other generalizations 
of the Bose-Hubbard model are of also possible by incoporating the influence of disorder. 
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